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ABSTRACT 

We study the thermal structure and evolution of magnetars as cooling neutron stars 
with a phenomenological heat source in an internal layer. We focus on the effect of 
magnetized {B > 10^** G) non-accreted and accreted outermost envelopes composed 
of different elements, from iron to hydrogen or helium. We discuss a combined effect 
of thermal conduction and neutrino emission in the outer neutron star crust and 
calculate the cooling of magnetars with a dipole magnetic field for various locations 
of the heat layer, heat rates and magnetic field strengths. Combined effects of strong 
magnetic fields and light-element composition simplify the interpretation of magnetars 
in our model: these effects allow one to interpret observations assuming less extreme 
(therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in 
their cores, can have higher thermal surface luminosity. 

Key words: dense matter — stars: magnetic fields — stars: neutron - neutrinos. 



1 INTRODUCTION 



We continue theoretical studies (iKaminker et al. l2006bl . 
hereafter Paper I; also see IKaminker et al.ll2006al . l2007l ) of 
persistent thermal activity of magnetar candidates - com- 
pact X-ray sources which include soft gamma repeaters 
(SGRs) and anomalous X-ray pulsars (AXPs). The magne- 
tars are thought to be warm, isolated, slowly rotating neu- 
tron stars of age t < 10 ^ yr with super strong magnetic fields 
B > 10" G (see, e.g.. lWoods fc Thompson 2.006, for a re- 
view) . Following many a uthors (e.g. , Colpi. Geppert fc Pagel 
I2OOOI : iThompsonI I2OOII : IPons et all 120071 '), we assume that 
the high level of magnetar X-ray emission is supported by 
the release of the magnetic energy in their interiors. Al- 
though this assum ption is widespread, there are alterna- 
tive models (e.g., Chatteriee. Hernauist fc NaravanI I2OO0I: 
Alpar '2001'; lErtan et all l2007l: iThompson fc Beloborodovl 
2005; Bclobo rodov fc Thompsonll2007l '). 



In Paper I we studied the thermal evolution of mag- 
netars as cooling isolated neutron stars with a phenomeno- 
logical heat source in a spherical internal layer. We ana- 
lyzed the location and power of the source and compared 
our calculations with observations of SGRs and AXPs. We 
showed that the heat source should be located at densi- 
ties p < 4 X 10^^ g cm~^, and the heating rate should be 
~ 10^'^ erg s~^ to be consistent with the observational data 
and with the energy budget of isolated neutron stars. A 
deeper location of the heat source would be extremely in- 
efficient to power the surface photon emission, because the 
heat would be carried away by neutrinos. 



Here, we refine the model of the magnetar heat- 
blanketing envelope. In this envelope, the effects of super- 
strong magnetic fields are especially important. We analyze 
the blanketing envelope consisting not only of iron (Fe) but 
also of light elements (H or He) which can be provided by 
accretion at the early stage of magnetar evolution. Chem- 
ical composition and strong magnetic fields do affect ther- 
mal conduction in the blanketing envelope and the thermal 
structure of magne tars. In addition to th e re sults of Paper I 
(and those of Potekhin fc Yakovlevl200ll and lPotekhin et"al] 
I2003D . we take into account neutrino energy losses in the 
outer crust of the neutron star, which can also be important 
(|Potekhin. Chabrier fc YakovlevllioOTi ). 



2 OBSERVATIONS 

For the observational basis, we take seven sources: three 
SGRs and four AXPs listed in Table [1] We present their 
ages t, efi^ective surface temperatures (redshifted for a 
distant observer), and redshifted thermal luminosities L^. 
The estimates of spin-down ages t, blackbody effective sur- 
face temperatures , and non-absorbed thermal fluxes are 
taken from the SGR/AXP online Catalog maintained by the 



McGiU Pulsar Gro up H All references in Table [J except for 
Rho fc Petr3l997l and rTiengo. Esposito fc MereghettilboOSl . 



are taken from that Catalog. We do not include SGR 1627- 
41, the faint X-ray pulsar CXO J164710. 2-455216 (e.g.. 



http:/ /www. physics. mcgill.ca/ ~pulsar/magnetar/main.html 
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Figure 1. Observational magnetar box (the shaded rectangle) 
of blackbody luminosity limits, (see text and Table [TJ, of 

seven magnetars, supplemented by observational limits on 
for eleven isolated neutron stars (marked by asterisks). The data 
are compared with theoretical cooling curves L^{t) for a 1.1 Mq 
neutron star with the dipole magnetic field Bp = 5 X 10^'' G 
(at the magnetic poles) and no internal heating, either without 
superfluidity (the solid line) or with strong proton superfluidity 
in the core (the dashed line SF). 



iMuno et al.ll200^ ) a nd the unconfirm ed AXP candidate AX 
J1845. 0-0258 (e.g.. iTam et aLlbOOel ). because their ages t 
are unknown (as well as the thermal luminosity of 
SGR 1627-41). We also do not include two AXPs, XTE 
J1810-197 a nd IE 1048.1-5937. Their h igh pulsed fractions 
(> 40%, e.g.. lWoods fc Thompsonll2006l) indicate that their 
flux comes from a small fraction of the surface, incompati- 
ble with the model considered here. In addition, we have ex- 
cluded the AXP 4U 0142-f 61 whose observed thermal emis- 
sion can be attributed to a circ umstellar dusty disk, rath er 
than to a neutron star surface (iDurant fc Kerkwiikll200^ . 

In Fig. [T]we plot the blackbody surface luminosity 
of the selected sources versus their ages. The current data 
are uncertain and our cooling models are too simplified to 
explain every source by its own cooling model. Instead, we 
will interpret magnetars as cooling neutron stars belong- 
ing to the "magnetar box," the shaded rectangle in Fig. [T] 
(which reflects an average persistent thermal emission from 
magnetars, excluding bursting states). 

The thermal luminosity limits can be obtained as 



4-kD fthxcex- 



(1) 



Here, /thx is a non-absorbed thermal flux detected from 
a source (in a certain E1-E2 X-ray energy band), D is a 
distance to the source and 

-1 



Qx = 



15 



da; 



exp(x) 



(2) 



is the bolometric correction, with xi^2 ~ -E'i,2/fcBT. In par- 
ticular, for the 2-10 keV band, we have ax ~ 2-4. Ther- 
mal fluxes fthx should be inferred from observations with 



account for the fractions BB/(PL+BB) of the blackbody 
{BB) components in the appropriate power law plus black- 
body {PL+BB) spectral flts (references are listed in Table 
[1]). The same fits provide the effective surface temperatures 
and the apparent radii Re's of emitting regions. The 
radii are defined in such a way that 



(3) 



where a is the Stefan-Boltzmann constant. The radii R'bb 
are typically smaller than the expected neutron star radii 
indicating that thermal emission can originate from some 
fraction of a neutron star surface. For instance, R^g «2.4 
and 5.5 km, for SGR 1806-20 and IE 1841-045, respectively. 
Although the actual surface temperature may strongly vary 
within the emission region, spectral fits give a single (surface 
averaged) value. The luminosities in Table [T] and 
Fig. [1] are mainly obtained from Eq. (|3]) using the values of 
Ts°° and R'bb presented in cited papers. 

CXOU J010043. 1-721134 is the only source from our 
collection, whose cumulative spectru m cannot be fit wit h 
a power-low plus blackbody model (|Tiengo et all I2OO8I ). 
iTiengo et all l|2008l ) fitted it by a sum of two blackbody 
components. The blackbody temperature of the softer com- 
ponent is given in Table [T] The radius of the corresponding 
emission region, R'bb = 12. ll^ ^ km, is comparable with the 
theoretical neutron star radius, although the pulsed fraction 
in the 0.2-6 keV energy range is rather high, 32 ± 3%. The 
harder component corresponds to — 7.89 io'gi MK and 
R'b'bs ~ ^-7 -o's km, meaning probably a hot spot on the 
neutron star surface. We define the thermal luminosity 
of this source as a sum of thermal luminosities of both black- 
body components. 

Radiation from five of the seven selected sources has the 
overall pulsed fraction < 20 %; the pulsed fr action for three 
of them is < 10% (e.g.. Woods fc Thompson 2006). This in- 
dicates that the thermal radiation can be emitted from a 
substantial part of the surface (although the pulsed fraction 
is lowered by the grav itational bending of light rays; e.g., 
iPavlov fc ZavlirJ I2OO O' and references therein). Two other 
magnetars from TablefT] CXOU J010043. 1-721 134 and AXP 
IRXS J170849.0-40091 0, have pulse fr action < 40% in the 
0.5-2.0 keV band (e.g.. lRea et al.ll2005l ). 

The majority of magnetar ages listed in Table [1] are 
characteristic spin down ages. For three sources, we adopt 
the ages of their host supernova r emnants (SNRs): t ^ 5 kyr 
for SGR 0526-66 in SN R N49 (|Kulkarni et al.l |2003| . also 
see lVancur a et al."l992^: t ~ 2 kyr for IE 1841-45 in SNR 
Kes 73 IVasisht fc Gotthclf 199 7.); and t 19 kyr for IE 
2259+586 in SNR CTB 109 (|Rho fc Petrelll997l ). To specify 
the left and right boundaries of the magnetar box in Fig. [1] 
we introduce, somewhat arbitrarily, the uncertainties by a 
factor 2 into the ages t. 

In our previous work we have compared simulations of 
magnetar cooling with the data on the effective surface tem- 
perature TI°°. Here, in contrast, we use the data on L^, 
which seem more robust. Note that explaining the data ei- 
ther on Ts°° or on with our cooling models is not entirely 
self-consistent. According to observations of all sources, but 
CXOU J010043.1-721134, thermal emission originates from 
some fraction of the magnetar surface while our cooling mod- 
els give thermal radiation from a large fraction of the surface. 
If we regarded (like in Paper I) the temperatures given 
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Table 1. Observational limits on ages t, effective surface temperatures Tj,°°, and blackbody luminosities LJ° of magnetars 



N 


Source 


t 




lgLJ° 


Refs.=) 






kyr 


MK 


erg s"-"- 




1 


SGR 1806-20 


0.22 


7.5+0.8 


35.1 +0.4 


MOO, M05 


2 


SGR 1900+14 


1.1 


5.0+0.6 


34.9 +0.5 


WOl 


3 


IE 1841-045 


2.0 


5.1+0.2 


35.15 +0.15 


VG97, M03 


4 


SGR 0526-66 


5.0 


6.2+0.7 


34.8 +0.4 


K03 


5 


CXOU J010043. 1-721134 


6.8 


3.5+0.2 


35.4 +0.2 


T08 


6 


IRXS J170849.0-400910 


9.0 


5.3+0.1 


34.6 +0.1 


R05 


7 


IE 2259+586 


19 


4.77+0.05 


34.4 +0.2 


RP97, W04 



°) Ages of SNRs (see text for references) 

The soft component of the double blackbody (BB) spectral model (see text for details) 
'^^ Total BB luminosity of both components of the double spectral model 

From Eq. with D = 3.0 + .5 kpc, the flux /thx in the 2-1 keV en ergy band, arid Qx ~ 3.8 
'=) MOO - Mcr eghetti et al.l ll2000h: M05 -iMereghetti et al.l (|2005^: WOl -IWoods et al.l ll200lh: yG97 - IVasisht fc Gotthelj lll997l'); M 03 
- iMorii et al. ( 200 ; KO3'" - iKulkarni et al.l | |2003^ : T08 - iTiengo et al.l 1 I2OO8I) ; R05 - iRea et ahl ll2005l) : RP97 - iRho fc Petrd l ll99it) : 
W04 - IWoods et all l|2004h 



by spectral fits (see Table [TJ as surface- averaged effective 
temperatures and calculated using Eq. ([Sjl witli values 
of Re's realistic for neutron stars, we would obtain notice- 
ably larger tiian those provided by the observations (the 
magnetar box would raise in Fig.[l|. On the contrary, match- 
ing the theory with the data on (as in the present pa- 
per) gives lower surface-averaged temperatures than the 
temperatures inferred from spectral fits. 

We expect that the theory will be improved in the fu- 
ture by constructing more advanced models of magnetars 
- for instance, with highly nonuniform sources of internal 
energy release. On the other hand, current interpretation of 
magnetar observations is far from being perfect. It would 
be a challenge to construct new models of thermal radia- 
tion from strongly magnetized neutron stars and use them 
(rather than blackbody models) to interpret the data. In 
this case, by analogy with employing hydrogen atmosphere 
models for describing thermal radiation from ordinary neu- 
tron stars, we expect to get higher Rbb (closer to the real 
neutron star radius) and lower . Moreover, thermal radi- 
ation emitted from a magnetar surface can be strongly dis- 
torted by magnetosp heric effects (e.g., iLvutikov fc Gavriill 
|2006I :Ir ea et al ] |2008D . This can greatly complicate the prob- 
lem of inferring and L|° from the data. 

Fig. [T] shows also observational limits for eleven 
ordinary iso lated neutron st a rs. Th e data are taken from 
Table 1 of iKaminker et al. 1 2006ah with a f ew changes 

20081 ). Following ISlane et al.1 



described by lYakovlev et al] 



2008*) and'Shiba nov et al.l l|2008t ). we have enlarged the age 
range of the pulsar J0205+6449 (in the SNR 3C 58) to its 
characteristic age of 5.4 kyr. We have excluded one young 
and warm source, IE 1207.4-5209, because of the prob- 
lems of interpre tation of its spectrum. We use the results of 
iHo et al.1 (|2003) for the neutron star RX J1856. 5-3754. The 
authors employed the magnetic hydrogen atmosphere model 
and obtained T^°° = (4.34 ± 0.02) x lO'' K and the appar- 
ent neutron star radius « 17 km (at the 68% confidence 
level for the fixed distance D = 140 pc). Taking into account 
a large scatter of distance estimates for RX J1856. 5-3754 



(|Walter fc Lattimej |2002| : Ivan Kerkwiik fc Kaplan! bOOTl ). 

we have added 10% error bars to the latter values of . 

Fig.[T]shows two typical cooling curves L^{t) for a low- 
mass neutron star (M = 1.1 A^q) without internal heating 
and with a strong dipole magnetic field {Bp = 5 x 10^"* G 
at the magnetic poles). The solid line is for a non-superfluid 
star, while the dashed line SF assumes a strong proton su- 
perfiuidity in the stellar core. This superfluidity suppresses 
neutrino emission in the core a nd thereby increases a t 
the neutrino cooling stage (e.g.. lYakovlev fc Pethic^l2004l ). 

Let us stress that the surface temperature of these stars 
is highly nonuniform; the magnetic poles are much hotter 
than the equator. In all figures we plot the total bolometric 
luminosity prod uced by the fiux inte grated over the stel- 
lar surface (e.g.. IPotekhin et al.ll2003l ). The observations of 
ordinary neutron stars can be ex plained by the cooling the- 
ory without any reh eating (e.g., lYakovlev fc Pethickl 12004 : 
lYakovlev et al ]|2008l ). The magnetars are much hotter (more 
luminous) than the ordinary cooling neutron stars; their ob- 
servations imply that they have additional heat sources. As 
in Paper I, we assume that these sources are located inside 
magnetars. 



3 PHYSICS INPUT 

We have performed calculations w ith our cooling code 
(|Gnedin. Yakovlev fc PotekhinlliooH ). which simulates the 
thermal evolution of an initially hot star via neutrino emis- 
sion from the entire stellar body and via heat conduc- 
tion to the surface and thermal photon emission from the 
surface. To facilitate calculations, the star is divided into 
the bulk i nterior and a thin outer heat-blanketin g enve- 
lope (e.g., iGudmundsson. Pethick fc EpsteinI Il983l ) which 
extends from the surface to the layer of the density p — 
Pb ~ 10^" — 10^^ g cm~'^; its thickness is a few hundred 
meters. 

In the bulk interior (p > Pb), the code solves the full set 
of thermal evolution equations in the spherically symmetric 
approximation. The standard version of the code neglects 
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4 X 

10^" g cm"^ 



the effects of magnetic fields on thermal conduction and 
neutrino emission. In the present version, we have included 
neutrino-pair electron synchrotron radiation in a magnetic 
field B, that was neglected in Paper I. 

In th e blanketing envelope , the updated version of the 
code (see IPotekhin et al ] |2007l . for details) uses a solution 
of stationary one-dimensional equations of hydrostatic equi- 
librium and thermal structure with radial heat transport, 
aniso tropic temperature distribu tion, and a dipole magnetic 
field (jGinzburg fc Ozernoi|[l963 ). It takes into account neu- 
trino emission and possible heat sources in the envelope. 
The solution, applied to different parts of the envelope with 
locally constant magnetic fields, yields temperature profiles 
slowly varying from one part to another. For a given T — Th 
at /9 = pb, we calculate the thermal flux emergent from 
different parts of the surface. Integrating it over the sur- 
face, we obtain the total photon luminosity Lg = A-koR^T^, 
where Ts is the effective temperature properly averaged over 
the stellar surface and R the circumferential neutron star 
radius. Redshifting then for a distant observer, we have 
Tr = Te s/l-r^/R and = (1-rg/i?) (rg = 2GM/c^ 
being the Schwarzschild radius). 

In the present calculations, we set either pb 
10^^ g cm~^ or (in the majority of cases) pb 
(see Sect. |4]). We consider blanketing envelopes composed 
of ground-state or accreted matter. For the ground state 
matter, the compositi on is iron up to p = 10^ g cm~^ and 
heavier elements (e.g., iHaensel. Potekhin fc Yakovlevll2007l ) 
at higher p (it will be called the Fe composition). As an 
alternative, we have studied fully accreted envelopes com- 
posed successively of H, He, C, O up to maximum p and 
T, where these elements can survive against pycno- or ther- 
monuclear burning, and then the Fe composi tion (we use the 
same structure of the accreted envelope as in lPotekhin et af] 
We have also considered accreted envelope models 
with H replaced by He, but found no significant difference. 

An anisotropy of heat conduction, produced by strong 
magnetic fields, essentially modifies the temperature dis- 
tribution in the blanketing envelope (which is included 
in our calculations). The anisotropy of heat conduc- 
tion can also create anisotropic temperature distribu- 
tion at p > Pb, particularly, in the inner crust (that 
is not included in our code). Such situations should 
be simulated with a two-dimensional coo l ing c o de, a s 
has been done by ICeppert. Kiiker fc Pagel (|2004 |2006| ): 
iPerez-Azorm. Mirrales fc Pond ||2006|) ffor stationary cas es) 
and most recen t ly bylAguil era. Pons fc Miralle^ ('2008a''P); 
lAguilera et"al] l|2009l ): iPons, Mirallcs fc Gcppcrt (2009). 
The effect strongly depends on the values of thermal con- 
ductivity across the magnetic field. If one restricts oneself 
to the thermal conductivity of strongly degenerate elec- 
trons, the anisotropy of heat conduction (the ratio of ther- 
mal conductivities along and across the field) in a deep 
and not very hot magnetized crust can be huge (because 
of the strong magnetization of electrons which are mainly 
moving along magnetic field lines). However, in a hot crust 
(with increasing temperature) the electron magnetization 
and heat conduction anisotropy weaken. Moreover, the con- 
ductivity of phonons (lattice vibrations of Coulo mb crys- 
tals of atomic nuclei; IChugunov fc Haensell |2007| ) and vi- 
brations of superflu id neutron liquid (superfiuid phonons; 
lAguilera et al.ll2009l ) in the inner crust can be much larger 



Table 2. Five positions of the heat layer, and the heat power 
W°° for the 1.4 Mq star with _ffo = 3 x 10^° erg cm-^ s"^ and 
t = 1 kyr 



No. 



pi (g cm 3) p2 (g cm ^) W°° (erg s ^) 



I 

II 
III 

IV 
V 



3 X 10^'^ 

10^2 
3 X 1013 
2 X 10^0 

10" 



10" 
3 X 10^2 

1014 
6 X 10"' 
3 X 10" 



4.0 X 10=^^ 
1.9 X 10^^ 

1.1 X 10^8 
3.3 X 10=*^ 
4.3 X 10^^ 



than the electron thermal conductivity across the field lines, 
washing out the temperature anisotropy and producing a 
nearly isotropic temperature distribution in the deep crust 
(although the electron conduction still dominates along field 
li nes). This tendency is quit e visible in recent calculations 
of lAguilera et al.1 (|2009l ') and lPons et all l|2009l ) and justifies 
our approach. 

At Tb > 10^ K the neutrino emission affects the ther- 
mal structure of the star. We calculate the neutrino emission 
in the crust (in the blanketing envelope and deeper) taking 
into account electron-positron pair annihilation, plasmon de- 
cay, neutrino bremsstrahlung in collisions of electrons with 
atomic nuclei, and synchrotron radiation of neutrino pairs by 
electrons (e.g., Yakovlov ct al. 20013- The pair annihilation 
is relatively unimportant and can be neglected. Very strong 
magnetic fields in the outer crust of magnetars can modify 
the plasmon decay and bremsstrahlung neutrino processes. 
Such modifications have not been studied in detail; they can 
be important for the magnetar physics but are neglected 
here. Evidently, the synchrotron emission does depend on 
the magnetic field, which we take int o account explicitly; it 
is important at B ~ IQi^'-lOi^ G (cf. iPotekhin et al.ll2007l '). 

Following Paper I, we introduce an internal phenomeno- 
logical heat source located in a spherical layer, pi < p < p2. 
The heat rate H [erg cm~'^ s~^] is taken in the form 

H = Hoe{p)eM-t/T), (4) 

where Ho is the maximum heat intensity, O ~ 1 in a density 
interval pi < p < p2, and Q vanishes outside this interval 
(at p <C pi or p ^ P2), f is the star's age, and r is the 
e-folding decay time of the heat release. An exact shape of 
0(p) is unimportant for Ts, provided that we fix the total 
heat power 



W°°{t) 



I 



dl/e'* H, 



(5) 



where dV is a proper volume element and $ is the metric 
function that describes gravitational redshift. An illustration 
of this independence is given in Fig. [S] which displays tem- 
perature distributions in the crust for two different shapes of 
0(p). We see that a change of a shape of 0(p) leaves almost 
intact the thermal structure at p ^ pi (and, therefore, it 
cannot affect thermal radiation). 

According to Paper I, only the values r ~ 10* — 10^ yr 
can be consistent with the magnetar box (Sect. [2]). Longer r 
would require too much energy (Sect.[S|. In this paper, we 
take r = 5 X lO"* yr. 

We employ the same EOS in the neutron star core 
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Figure 2. Temperature profiles along the magnetic axis of the 
neutron star with M = 1.4 M©, B = 10^^ G, Ho = 3 X 10^" 
erg cm~^ s~^, and two forms of the heat source profile, shown in 
the inset: smooth (solid lines) and piccewise-constant (dot-dashed 
lines). 



as in Paper I. I t is th e model denoted as APR III 
by iGusakov et al.1 (120051): it is based on the EOS of 
lAkmal. Pandharipande fc Ravenhalll l|l998l ). According to 
this EOS, the core consist of nucleons, electrons, and 
muons. The maximum neutron star mass is M = 1.929 Mq. 
The powerful direct Urea process of neutrino emission 
l|Lattimer et al.l Il99ll ) is allowed only in the central ker- 
nels of neutron stars with M > 1.685 Mq (at densities 
p > 1.275 X 10^^ g cm"'*). 

We use two neutron star models, with M = 1.4 Mq and 
1.9 Mq. The former is an example of a star with the standard 
(not too strong) neutrino emission in the core (the modi- 
fied Urea process in a non-superfluid star). In this ease R = 
12.27 km and the central density is pc = 9.280x 10^* g cm"^. 
The latter model {R = 10.95 km, pc = 2.050 x 10^^ g cm"^) 
gives an example of a star whose neutrino emission is en- 
hanced by the direct Urea process in the inner core. 

Five examples of heat layer locations, pi, p2, are given 
in Table[2l Three of them (I, II, and III) were used in Paper I. 
Let us remind that the outer crust has a thickness of a few 
hundred meters and a mass of ~ 10"^ Mq; the inner crust 
can be as thick as ~ 1 km and its mass is ~ 10~^ Mq, 
while the core has radius ~ 10 km and contains ~ 99% 
of the stellar mass. All five heat layers are relatively thin. 
The layers I, IV, and V are located in the outer crust; the 
layers II and III are at the top and bottom of the inner 
crust, respectively. For illustration, in Table [2] we present 
also the heat power W°° calculated from Eq. ((SJ for the 
five layers in the 1.4 Mq star of age t — 1000 yr at Ho = 
3 X 10^" erg em"^ s"^ 



4 THERMAL STRUCTURE OF 

HEAT-BLANKETING ENVELOPES 

Fig. |3] shows the temperature and radial thermal fiux pro- 
files in local parts of the heat blanketing Fe envelope (p ^ 
Pb = 4 X 10" g cm"^, Tb = lO'' K) of the 1.4 Mq star with 
a locally uniform magnetic field B — 10^"' G directed at two 
angles to the surface normal, 9b = and 9b = 75°. One can 



compare the thermal structure of the blan keting envelope 
without he ating (short-dashed lines; e.g., IPotekhin et al.l 
|2003| . |2007| ) and with the heat source (solid lines) of the 
intensity Ho = 3 x 10^*^ erg em~^ s~^, located in layer I 
(Table il). 

Long-dashed lines in Fig. |3] are obtained by solving 
the full se t of one-dimens i onal equations for the blanketing 
envelope (|Potekhin et al.l l2007 t) including the heat source 
but with the magnetic field artificially switched off at p > 
10^" g em~'^. These calculations simulate the model used 
in Paper I, where a strongly magnetized blanketing enve- 
lope with Pb = g em~^ was matched to the interior, in 
which the magnetic field effects were ignored. We see that at 
different 6b the long-dashed r(p) curves only slightly differ 
from the solid ones. Such a difference is more pronounced 
near the heat layer (the bottom left panel) but becomes in- 
visible at lower p. We have obtained a significant difference 
only in a narrow range of field directions 9b ~ 90°. However, 
in these eases the one-dimensional (radial) model becomes a 
poor approximation. The curvature of magnetic field lines in 
a more realistic model should increase the surface tempera- 
ture in a narrow equatorial zone of width < V Rh (along the 
surface) , where h is the thick ness of the heat-blanketing en- 
velope (jPotekhin et al.ll2007h . An additional increase of the 
surface temperature in the above region can be p rovided by 
ion heat conduction (|Chugunov fc Haensell l2007l) . The tem- 
perature raise will reduce the indicated difference between 
the solid and long-dashed lines. 

We have obtained similar results for local radial heat 
flux F{p) shown in the bottom-right panel of Fig. [S] The 
flux changes its sign inside the heat layer. The flux at 
P 10^^ g em~^ flows into the stellar interior, where the 
heat is radiated away by neutrinos (Paper I). The solid 
and long-dashed curves F{p) are also indistinguishable at 
p < 10^" g em~'^. Calculations show that the convergence of 
two types of the curves is violated only at B > 10^® G. 

We have verified that our calculations for the blanketing 
envelope with pb = 4 x 10^^ g cm^'^ (including the heat 
layer) properly match those with pb = lO'^" g em~^ (with 
the same heat layer being outside the blanketing envelope). 
These results justify the choice of the blanketing envelope 
with pb = 10^° g em~^ in our further calculations (at least 
with B < 10^^ G). 

Fig. 2] shows the average surface temperature Ts of the 
Af — 1.4Mq star versus Tb. We assume a dipole mag- 
netic field with Bp = 10^^ 10^^ and 10^^ at the magnetic 
poles. We consider our Fe and accreted envelopes as well 
as an accreted envelope with all hydrogen replaced by he- 
lium. One can see an appreciable increase of Ts with the 
growth of Bp above 10^* G because of the cumulative ther- 
mal conductivity enhancement. The most significant effect 
of the accreted envelopes is a systematic increase of Ts at 



been studied earlier ( 


Potekhin. Chabrier fc YakovlevI Il997l: 


iPotekhin & YakovlevI 


20011: IPotekhin et al.ll2003l). Here, we 



have verified that it is important for all magnetic fields Bp 
of our interest; it results from the thermal conductivity en- 
hancement in large areas of the accreted envelope near the 
magnetic poles. More details on the Tb — Ts relation for Fe 
envelopes are given in the Appendix. 

The right panel in Fig. |4] shows that replacing hydro- 
gen by helium in the accreted envelope at Tb > 10^ K does 
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Figure 3. Temperature T(p) and radial thermal flux F(p) profiles in local parts of the outer crust of the 1.4 Mq neutron star (^®Fe 
envelope; = 4 X lO'^^ 

the surface normal, 9g = (thin lines) or dg 
strips; the heat intensity is = 3 x 10^" erg cm 



cm~'^ and fixed T^, = 10 K) with a locally uniform magnetic field B = 10 G directed at two angles 9g to 
75° (thick lines). The position of the heat layer (I in Table |2j is indicated by the shaded 
^ s~^; solid lines - accurate calculations, long-dashed lines - magnetic field is off at 



p > 10 " g cm~^; short-dashed lines — no heating. Ho = 0. Top: entire temperature profiles at p ^ Pb, the filled dots show the radiative 
surface (optical depth equals 2/3). Bottom: T(p) in units of 10^ K (left) and F{p) in units of 10^'^ erg cm~^ s~^ (right) in the vicinity 
of the heat layer. 




Figure 4. Surface temperature Tg, averaged over the surface of the M = 1.4Mq neutron star, as a function of Tt, for the blanketing 
envelope with p^, = 10^" g cm~^. Three families of curves 12, 15, and 16 correspond to the dipole magnetic fields with Bp = 10^^, 
10^^, and 10^® G, respectively. Left: Solid and dashed lines refer to the Fe and accreted envelopes; dotted lines are for the Fe envelopes 
neglecting neutrino emission. Right: Dashed lines - standard accreted envelopes; dash-and-dot lines 'acc (He)' - the accreted envelopes 
with hydrogen replaced by helium; dotted lines - standard accreted envelopes neglecting neutrino emission. 
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Figure 5. Average surface temperature Ts as a function of Tj, 
for tlie 1.4Mq (thick lines) and 1.9Mq (tliin lines) stars; solid 
and dashed lines refer to the same envelopes as in Fig. |4] labels 
12 and 16 refer to the field strengths Bp = lO^^ and W^^ G. 

not affect T^. The insensitivity of Tg to t his replacement 
for n on-magnetized neutron stars is known (jPotekhin et al.l 
1 19971 ). Here, we have checked this property for magnetars. 
Only at Bp ^ 10^* G and relatively low temperatures 
Tb < 10® K the surface temperature Ts for the He enve- 
lope goes slightly higher. Generally, such an effect is un- 
usual for ions with higher which are better scatterer of 
electrons (cf. IPotekhin et~allll997l . l2003l ). However, it oc- 
curs in the superstrong field because of deeper localization 
of the radiative surface for the He envelope, as a result of 
a lower energy plasma-frequ ency cutoff in the Rosseland 
opacity jPotekhin et al.ll2003l ) for smaller Z/A. Anyway at 
Tb ^ 10* K replacing hydrogen by helium does not affect 
the thermal insulation of the blanketing envelope. 

Also, Fig. |4]shows the effects of neutrino emission in the 
blanketing envelopes of different composition at Tb > 10^ K. 
One can see that the neutrino emission limits the growth of 
T; with increasing Tb (cf. [Potckhin et al. 2007). 

Fig. [5] shows the dependence of Ts on Tb for the 1.4 Mq 
and 1.9 Mq stars with Fe and accreted envelopes. In a wide 
range of magnetic fields, we obtain systematically higher 
surface temperatures Ts (at the same Tb) for the massive 
star as a result of smaller radius (see Sect. [3} and thinner 
blanketing envelope (A_Rb = R — Rh ~ 170 m for AI — 
1.4 Mq versus AJ?b « 70 m for M = I.9M0). At a given Tb, 
the surface temperature scales approximately as Ts oc jr^''*, 
where (j ~ GM / (B? \/l — r^^/R) is the surface gravity (e.g., 
IVentura fc Potekhirj|200ll , and references therein). 



5 COMPARISON WITH OBSERVATIONS 

A warm cooling neutron star with a powerful internal heat- 
ing (Sect. [3]) quickly (in t < 10 yr) reaches a quasi-stationary 
state regulated by the heat source. The energy is mainly 



carried away by neutrinos, but some fraction is transported 
by thermal conduction to the surface and radiated away by 
photons; the stellar interior stays highly non-isothermal. A 
thermal state of the heat source and outer layers is almost 
independent of the physics of deeper layers. This means the 
thermal decoupling of the heat source and outer layers from 
the deeper layers. 

In general (Paper I) heating a warm neutron star from 
the core or the inner crust is inefficient for raising Ts. With 
the increase of the heat power, Ts saturates because of the 
neutrino emission. Similar saturation of Tg b y the neutrino 
emission was analyzed by IVan Ripej |im3) 

who studied 

thermal response of a neutron star to a steady-state heat- 
ing. In Paper I we concluded that the heat source should be 
located in the outer crust in order to heat the surface and 
be consistent with the neutron star energy budget. 

Fig. |6] is similar to Fig. 2 of Paper I. It is a reference 
figure for subsequent Figs.[7H9l The left panel of Fig.|6]shows 
the temperature profiles Ti(p) inside the 1.4 Af© star of age 
t = 1000 yr with the dipole magnetic field (Bp = 5x lO" G). 
Here, Ti(p) = T{p)e^ is the internal temperature redshifted 
for a distant observer, while T(p) is the local temperature 
at a given p. It is Ti that is constant throughout thermally 
relaxed (isothermal) regions of the star in General Relativity. 
The s ame temperature Ti(p) has been plotted in Paper I 
and in lKaminker et al.l (|2007| ) (denoted there as T{p)). We 
consider three locations of the heat layer (I, II, and HI in 
Table and two intensities, i^o = 3 x 10^^ and 3 x 10^" 
erg cm~^ s~^. In all the cases the stellar core is colder than 
the crust, because of intense neutrino cooling in the core. 
Pushing the heat source deeper into the crust we obtain a 
colder surface because of more efficient neutrino cooling. 

The right panel of Fig. [S] shows cooling curves. Nearly 
horizontal parts of the curves at < r = 5 x 10* yr (Sect. 
[3| and their later sharp drops confirm that the surface ther- 
mal luminosity is solely maintained by internal heating. The 
initial parts {t < 100 yr) of two cooling curves for the heat 
layers II and III (at Hq ~ 3 x 10^^ erg cm~^ s~^) demon- 
strate the end of the relaxation to quasi-stationary thermal 
states. We can reconcile the cooling curves at lower heat in- 
tensity with more luminous sources in the magnetar box by 
placing the heat source in the outer crust. To minimize the 
energy consumption (see below) we employ the outer heat 
layer I in Figs. [T] |S]and the layers I, IV, and V in Fig. (5] 

The left panel in Fig. [7] shows the temperature profiles 
Ti(p) inside the 1.4 Mq star of the age t = 1000 yr with the 
heat source in layer I calculated for both the Fe and accreted 
blanketing envelopes. For comparison, we take three levels of 
the heat intensity, i/o = 3 x 10^® (thin lines), i/o = 3 x 10^^ 
(intermediate lines), and 3 x 10^" erg cm~'^ s~^ (thick lines), 
and two magnetic fields (Bp = 10^'' and lO^'' G). The right 
panel in Fig. [7]demonstrates the appropriate cooling curves. 

Under the heat-blanketing envelope (at p > pb), we 
take into account the magnetic field effects only by includ- 
ing the synchrotron neutrino radiation and (indirectly) the 
heat source that can be provided by magnetic fields. The 
synchrotron emissivity is calculated by putting B = Bp. 
The neutrino synchrotron process in superstrong magnetic 
fields (10^** — 10^^ G) lowers the temperature profiles in the 
stellar interior; this effect is more pronounced for stronger 
fields (cf. solid lines for Bp = 10^* and dot-dashed lines for 
Bp — 10^"' G in the left panel of Fig.[7ll. The accreted matter 
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Figure 6. Left: Temperature dependence Ti(p) in the 1.4 Mq neutron star of age t = 1000 yr with three different positions I, II, and III 
(Tabled of the heat layer (hatched rectangles) and two levels of the heat intensity Hq = 3 x 10^® erg cm"-^ s~^ and 3 X 10^*^ erg cm~'^ s~^ 
for the Fe blanketing envelope and the dipole magnetic field with Bp = 5 X 10^^ G at the poles. In the upper part we indicate the density 
regions appropriate to the outer crust, inner crust and the core of the star. Right: Cooling curves in comparison with the magnetar box. 
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Figure 7. Left: Same as in the left panel of Fig. [6] but for one heat layer I, three heat intensities, and two magnetic fields (Bp = 10^* 
— solid and long-dashed lines; and Bp = 10^^ G - dot-dashed and short-dashed lines). Solid and dot-dashed lines correspond to the 
Fe heat-blanketing envelope, dashed lines — to the accreted envelope; thick, intermediate, and thin lines are for the heat intensities 
Ho = 3x lO^s, 3 X 10i9, and 3 x 10^0 erg cm ^ s ^, respectively. Right: Cooling curves for these models. 



in the blanketing envelope also reduces the temperature at 
P ^ 10^^ g cm~^ because of higher heat transparency of the 
accreted envelope (higher heat flux to the surface) . 

Fig. [7] demonstrates that is mainly regulated by 
the blanketing envelope. The combined effect of a super- 
strong magnetic field and an accreted envelope apprecia- 
bly increases (also see Fig. |4)|. The stronger heat- 
ing {Ho = 3 X 10^° erg cm"^ s"^) at Bp = 10^^ G pro- 



duces too warm magnetar envelope for any composition and 
gives larger than required by the magnetar box. In a 
lower field, Bp — 10^'* G, the accreted envelope is also too 
warm but the Fe envelope is cooler and better consistent 
with the magnetar box. The weaker (intermediate) source 
{Ho — 3x 10^^ erg cm"'' s~^) also overheats the accreted en- 
velopes at both Bp = 10^^ G and Bp = 10^* G (intermediate 
short- and long-dashed lines); in this case the cooling curves 
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for Fe envelopes better match the data. However, only the 
weakest chosen heat intensity {Ho = 3 x 10^* erg cm~^ s~^) 
is capable to cover the lower part of the magnetar box. 

Fig. [S] shows even more pronounced effects of the mag- 
netar magnetic fields and accreted envelopes for the 1.9 M© 
star with two levels of heat intensities in the layer I. The 
left panel gives the temperature profiles Ti{p) a.t p > pb — 
10^" g cm~'^. They are noticeably lower than the correspond- 
ing profiles in the 1.4 Mq star (because of the direct Urea 
process that operates in the inner core of the massive star). 
However, in the case of Ho = 3 x 10^" erg cm~^ s~^ the 
effect of the magnetized accreted envelope overrides that of 
the rapid neutrino cooling of the massive core (owing to ther- 
mal decoupling of the surface from the core). Comparing the 
right panels of Figs. [7]and|8l we see that at Ho = 3 x 10^" 
erg cm~^ s~^ and t < 5 x 10* yr the thermal luminosity of 
the heavier 1.9 A/q star is higher than of the 1.4 star 
(because of higher Ts, see Fig. [5}, making the heavier star 
too hot. As discussed in Paper I, an intense heating in the 
outer crust of a massive star can outweigh fast neutrino cool- 
ing in the inner core. This effect is very unusual for ordinary 
cooling stars where massive stars are commonl y colder than 
low-mass ones (e.g.. lYakovlev fc Pethic^l2004l ) . 

In contrast, the effects of fast neutrino cooling in mag- 
netars are more essential at Ho = 3 x 10^® erg cm~'^ s~^. 
Tuning //o ~ (1 — 3) x 10^^ erg cm~^ s"^ and the chem- 
ical composition of the blanketing envelope, we can rea- 
sonably well explain the magnetar box. Note that the low- 
est heat intensity Ho = 3 x 10^^ erg cm""^ s~^ taken in 
Fig. [8] is ten times larger than the lowest heat intensity 
Ho = 3 X 10^* erg cm"^ s"^ taken in Fig. [71 but the corre- 
sponding cooling curves do not strongly differ. 

Earlier we (jKaminker et al.l l2007h have shown that 
strong variations of the thermal conductivity in the inner 
crust for the case of intense heating in layer I have no effect 
on the surface luminosity. Fig. [9] demonstrates the sensi- 
tivity of the thermal structure (the left panel) and cooling 
curves (the right panel) to artificial variations of the ther- 
mal conductivity in the outer crust. The left panel shows 
temperature profiles in the 1.4 Mq 1000 yr-old star with 
the dipole magnetic field (Bp = 5 x 10^'* G) for three lo- 
cations of the heat layer (cases I, IV, and V in Table [Sjl at 
Ho = 3 X 10^" erg cm~^ s~^; the right panel gives respec- 
tive cooling curves. The thick solid lines are the same as 
in Fig. [6] Other lines are calculated with the thermal con- 
ductivity modified in the density range 1.1 x 10^^ ^ p K 
4 X lO^'^ g cm""^ (the range is marked in the left panel of 
Fig. [9] by a dot-hatched rectangle). 

Thick dashed and dot-dashed lines are calculated with 
the thermal conductivity reduced by a factor of 10 (kcH — 
k/10). They illustrate a possible suppression of radial heat 
conduction (e.g., by a strong toroidal magnetic field in the 
outer crust). For the heat layers I (long-dashed lines) and 
IV (short-dashed lines) located closer to the bottom of the 
blanketing envelope, the conductivity reduction results in a 
sharper temperature drop inside the crust and in a cooler 
interior, with the tendency to the isothermal state. Taking 
the heat region V (thick dot-dashed lines), shifted to the 
inner edge of the layer with the reduced conductivity, we 
obtain qualitatively the same behavior of Ti(p), as in the 
cases II and III in Fig. [G] The thermal energy easier fiows 



inside the star than in the case I with normal conduction 
(solid lines). 

On the contrary, the enhanced thermal conductivity 
(thin dashed lines) produces a wide quasi-isothermal layer 
in the outer crust (the left panel of Fig. ^ and a photon 
surface luminosity (the right panel) that is nearly the same, 
as for the normal conductivity (the thick solid line) . In other 
words, L^{t) is slightly sensitive to a conductivity increase. 
Comparing the right panels of Figs. [9] and [7] (thick lines), we 
conclude that a conductivity increase at p > pb is incapable 
to rise , while an increase at p < pb can rise it. 

Finally, let us discuss briefly the energy budget of mag- 
netars. Following Paper I, we assume that the maximum 
energy of the internal heating is -Emax ~ 10^" erg (which 
is the magnetic energy of the star with _B ~ 3 x 10^® G 
in the core). Then the maximum persistent energy gener- 
ation rate is VFmax ~ S'max/T ~ 3 X 10^^ erg s~^. For ex- 
ample, let us take an 1.4 neutron star of age t <^ t 
with the heat source in layer I. For an intense heating with 
J/o ~ 3 X 10^" erg cm'^* s"^ we obtain W°° ~ VP'max (and, 
therefore. Ho cannot be larger). For a less intense heating 
with Ho ~ 3 X lO'^^ erg cm~"^ s~^ we have a more relaxed 
condition W°° ~ 0.1 Wmax (which would leave some energy 
for bursting activity of magnetars). 

It follows from Figs. [BHSl that the heating should be 
sufficiently intense to keep Lf on the magnetar values 
(~ lO'^^ erg s~^). However, for realistic magnetic fields 
B ~ (2 - 10) X 10" G, the maximum allowable heat 
intensity Ho ~ 3 x 10^" erg cm~"^ s~^ and accreted en- 
velopes, we have the thermal surface luminosity ~ 10'^'' 
erg s~^, noticeably higher than the luminosity of magne- 
tars (Fig. [?}• For Fe envelopes and the same heat inten- 
sity, we obtain ^ 3 x 10^^ erg s~^, consistent with the 
upper part of the magnetar box but giving the stringent 
energy budget {W°° ~ Wmax). Using a weaker heat inten- 
sity Ho ~ 3 X 10^^ erg cm~'^ s~^ and accreted envelopes, 
we obtain still greater thermal luminosity > 4 x 10"^^ 
erg s^^ , which is too high for the magnetar box but provides 
a reasonable energy budget. Finally, varying weaker heating 
rate Ho < 10^^ erg cm~^ s~^ and the chemical composition 
of the blanketing envelope (Fig. (Tjl we have the luminosity 
L'^ ~ 10^^ erg s~^, that is consistent with the magnetar 
box. 

Accordingly, the presence of accreted envelopes simpli- 
fies the explanation of magnetars as cooling neutron stars 
(in our model). Our results show that we can reconcile the 
theory with observations assuming the accreted envelopes 
and lower heat intensities. Ho ~ 10^® erg cm~^ s~^. Note 
that in all the cases the efficiency of heat conversion into 
the thermal radiation, L^/W^°°, is low but the accreted en- 
velopes increase it. For instance, assuming Ho ~ 3 x 10^^ 
erg cm~^ s~^ we have /W°° ~ 0.01 for Fe envelopes and 
LT/W^ ~ 0.1 for accreted ones. 



6 CONCLUSIONS 

We have analyzed the hypothesis that magnetars are 
isolated neutron stars with B > lO^** G, heated by a 
source localized in an internal layer. We have modelled 
the thermal evolution of magnetars, taking into account 
that their heat blanketing envelopes can be composed 
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Figure 8. Same as in Fig. [7]but for the 1.9 Mq star and two heat intensities Hq = 3 x 10^" (thin lines) and 3 X 10^" (thick lines). 
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Figure 9. The effects of thermal conduction in the outer crust on the thermal evolution of the 1.4 star with the same magnetic 
field as in Fig. |6]at Ho = 3 X lO'^" erg cm~'^ s~^. Left: Temperature profiles in the star at t =1000 yr. The hatched rectangles show the 
positions of the heat layers I, IV, and V, and the layer (labeled as K(,ff = ''^/lO), where the thermal conductivity is modified (see the 
text). Solid line is the same as in Fig.[6l thick short-dashed, long-dashed, and dot-dashed lines are for the thermal conductivity reduced 
by a factor of 10 and for the heat layers IV, 1, and V, respectively; thin dashed line is for the thermal conductivity enhanced by a factor 
of 10"^ and the heat layer I. Right: Cooling curves in comparison with the magnetar box. 



of light elements. Such envelopes can appear either 
due to a fallbac k accr et ion a fte r a supernova explosion 
(e.R. IChevaliej 1 19891 . 1 19961 : IChang. Arras fc Bildstenl 
I2OO4I ). probably with subsequent nuclear spal- 
lation reactions l|Bildsten. Salpeter fc WassermanI 



tion from a fossil disk (e.g., Chatteriee et al 


WanE. Chakrabartv & KaolarJ 


20061: Ertan et al. 


Romanova. Kulkarni fc Lovelace 


I2OO8I) or from the 



200c 



2007 



stellar medium (e.g.. [Nelson. Salpeter. fc WassermanI Il993l : 
iMorlcy 199i). 

Compared to Paper I, we have (i) included the effects 
of accreted envelopes, and (ii) changed the strategy of rec- 
onciling the theory with observations. We rely now on ob- 
servational limits of quasi-persistent thermal luminosities of 
magnetars; this lowers the magnetar box (Sect. [2]) and re- 
laxes theoretical constraints on the properties of internal 
heat sources. 
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The main conclusions are as follows: 

(1) The presence of light elements in the outer enve- 
lope of a magnetized neutron star can significantly increase 
the thermal conductivity and the thermal stellar luminosity 

(for a given temperature Tb at the bottom of the heat 
blanketing envelope). Similar conclusions have been made 
earlier for ordinary cooling neutron st ars w ith B < 10^'' G 
ijPotekhin et all 1 19971 : lYakovlev fc Pethickl booj . and refer- 
ences ther ein) as well as for s trongly magnetized cooling 
stars ("e.g. JPotekhin et al.llioO^ '). 

(2) The luminosity of the star with an accreted en- 
velope is insensitive to replacing all ac creted hydrogen by he - 
lium (as in ordinary neutron stars, see lPotekhin et ai]|l997l ). 
In particular, these results can be used for taking into ac- 
count rapid nuclear burning of hydrogen and accumulation 
of he lium in the outer part of the envelope (e.g.. lChang et al] 
|2004 ). 

(3) The combined effect of a superstrong magnetic field 
and an accreted envelope simplifies the interpretation of ob- 
servations of quasi-persistent thermal radiation from mag- 
netars using our model. We confirm the conclusion of Paper 
I that the our most favorable models require the heat source 
to be located in the outer crust (at p < 4 x 10^^ g cm~^). 
However, the presence of accreted envelopes allows us to take 
lower heat intensities Ho ~ 10^^ erg cm~'^ s~^ and place the 
heat layer slightly deeper in the stellar interior. 

(4) In accordance with Paper I, in all our successful 
models (with and without accreted envelopes), heating of 
the outer crust produces a strongly nonuniform tempera- 
ture distribution within the star. The temperature in the 
heat layer exceeds 10^ K, while the bottom of the crust and 
the stellar core remain much colder. The outer crust is ther- 
mally decoupled from the inner layers; thermal surface emis- 
sion is rather insensitive to the properties of the inner layers 
(such as the EOS, neutrino emission, thermal conductivity, 
superfiuidity of baryons). 

(5) The surface thermal luminosity is weakly affected 
by variations of the thermal conductivity in the outer crust 
below the heat blanketing envelope. Therefore, the effects 
of the magnetic field on the conductivity in the heat layer 
cannot greatly change the surface luminosity. The thermal 
surface radiation is mainly regulated by the heat source as 
well as by the magnetic field and chemical composition of 
the blanketing envelope. Nevertheless, our calculations can 
be improved by a more careful treatment of heat transport in 
the entire magnetized outer crust, at p > 10^' ^ g cm~'^, with 
different magnetic field con figurations (e.g., iGeppert et al.l 
l2006l : I Aguilera et al.|[2008ah . 

(6) Increasing the surface thermal emission of the 
star, which has a relatively high heat intensity {Ho ~ 
10'^° erg cm~'^ s~^) and an accreted envelope, is even more 
efficient if the star is massive (and possesses, therefore, thin- 
ner and more heat transparent crust). This effect is stronger 
than fast neutrino cooling due to direct Urea process that 
can be allowed in the core of a massive star. 

(7) The presence of an accreted envelope can raise the 
efficiency of heat conversion into the surface radiation. It 
can become as high as ~ 10% (compared to a maximum of 
~ 1% for Fe envelopes). This enables us to make our models 
more consistent with the total energy budget of heat sources 
in a neutron star. Now we can reduce the total energy to 



^tot ~ 10'** — 10*^ erg (instead of the previously assumed 
level of E^t ~ 10*^ - 10^° erg). 

Further observations as well as new models of magne- 
tar atmospheres are needed for more reliable interpretation 
of observations. The physics of internal heating is still not 
clear; it should be elaborated in the future. 
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APPENDIX A: Tb-Ts RELATION 

In Tables IA1I4A3I we present the relations between the 
temperature Tb at the inner boundary of the blanketing 
envelope (p = pb) and the non-redshifted effective surface 
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Table A2. Same as in Table IXT] but for Bp = lO^S G. 
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Table A3. Same as in Table lAl] but for = 4 X lO" g cm"-*. 
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Figure Al. The surface-averaged effective temperature Ts (left 
vertical axis, solid and dot-dashed lines) and the total heat flux 
Lrb through the inner boundary of the Fe blanketing envelope 
(right vertical axis, dashed and dotted lines) in the 1.4 Mq star 
for pb = 10^"^ g cm~^ (solid and dashed lines) and pb = 4 X lO'^^ 
g cm~'^ (dot-dashed and dotted lines), for B = 10^^ G (lower 
curve of each pair) and 10^^ G (upper curves). 

saturates at Tb > 10^ K if pb is placed at the bottom of 
the outer crust, pb = 4 x 10^^ g cm~^, but it does not 
saturate at lower pb = 10^" g cm~^. In the former case 
the saturation occurs because of the strong neutrino emis- 
sion at Pb > 10^" g cm~^ and high temperatures. We also 
plot the total outward heat flux through the boundary of 
the blanketing envelope, Lrb, in such a scale that corre- 
sponding curves match the Tb - Ts curves at low Tb. Then, 
the deviation of the Ts-curves from the Lrb-curves directly 
measures the total energy loss due to neutrino emission 
in the blanketing envelope (at p < pb). 



temperature Ts. In addition, we list the values of the out- 
ward radial heat flux Tb at p = pb. We assume no heat 
sources in the envelope and consider the Fe envelope in 
the 1.4 Mq neutron star of radius 7? = 12.27 km with 
the dipole magnetic field (Bp = 10^^ or 10^^ G at the 
magnetic pole). We present local values of Ts and Tb at 
magnetic poles, as well as surface-averaged values (av). 
Table \M\ refers to Bp = 10^^ G and pb = 10^° g cm"^; 
Table m to Bp = 10^^ G and pb = 10^" g cm"^; and 
Table HH to Bp = 10^^ G and pb = 4 x lO" g cm'^ 
Fig. lAll compares these Tb-Ts relations. We see that Ts 



